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Abstract 

"The plasma approximation" is, in the words of Pauli, "not even wrong," as it expresses a 
disbelief in the symmetries of the underlying gauge field theory. The electrostatic field in question 
is divergenceful yet has no sources or sinks, thus breaking Lorentz covariance explicitly, and has no 
place in a self-consistent theory. Utilizing it leads to inconsistencies in the equation of motion and 
prevents a proper, field-theoretic treatment of condensed matter in the plasma state. By exposing 
the fallacy of replacing one of the field equations with one of the Bianchi identities, the applicability 
of Gauss' law to plasma physics is assured. 
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All approximations are wrong; the question is, how wrong? "The plasma approximation" 
which is assumed applicable to low-frequency and steady-state phenomenon, popular among 
plasma scientists in various fields 3 an d tracing its roots beyond Braginskii 0, 4], is, in 
the words of Pauli, "not even wrong," as it expresses a disbelief in the symmetries of the 
underlying gauge field theory. (Many authors 5|,|6| use the term to mean what is better called 
"the plasma condition" [7], ie the number of electrons in a Debye sphere is large enough 
to effect charge shielding, a situation akin to the vacuum polarization found in quantum 

issue herein. What is at issue here is the presupposition 
of a diverging electrostatic field on scales much larger than the Debye length in the absence 
of an explicit supporting charge density.) The electrostatic field in question is divergenceful 
yet has no sources or sinks. Utilizing it leads to inconsistencies in the equation of motion 
and prevents a proper, field-theoretic treatment of condensed matter in the plasma state. 
As Griffiths says 9|, "like Gauss' law, Ampere's law [with Maxwell's correction] is always 
true.''' The reason is because these equations "may be regarded as a four- dimensional version 
of Poisson's equation." Lorentz covariance implies a symmetry in the inhomogeneous of 
Maxwell's equations which must be respected, as decoupling electrostatic fields from source 
charges is tantamount to decoupling magnetostatic fields from source currents. We reiterate 
the development of classical and quantum field theory 0, Q, H Q, Q and insist on its 
applicability to all media, including magnetized plasma. 

Let us first be pedantic with our terminology, as only formal definitions equate with formal 
mathematics. The term "quasineutral" relates to "the plasma condition" defined above, 
referring to the overwhelming tendency of electrons within a plasma to shift their positions in 
response to the electrostatic potential of the ion density so as to attenuate exponentially the 
Coulomb field, and is often taken as the definition of "the plasma approximation." However, 
the traditional term "neutral" already encompasses the implications of quasineutrality, as 



no medium (except that of neutron stars 14j|) remains neutral on scales sufficiently small 
to resolve isolated charges. Within neutral media, the microscopic field fluctuates wildly 
on the particle scale but averages out as the differential volume element grows 0, . In 
contrast, as Chen succinctly states [l|], "the plasma approximation is almost the same as 
the condition of quasineutrality discussed earlier but has a more exact meaning [and] is a 
mathematical shortcut that one can use even for wave motions." Specifically, "it is usually 
possible to assume n« = n e and V-E ^ at the same time." Indeed, most modern definitions 
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of the plasma magnetohydrodynamic equations 0, Q, 0, Q, Q, 

three-fourths of Maxwell's laws. Regardless of the name one gives the procedure, the idea 
that Gauss' law is inapplicable within a plasma runs contrary to well over a century's worth 
of experimental and theoretical analysis [23j and leads one into a host of difficulties in 
determining the physics of the plasma state. 



The 

Theory 



ol 



owing draws heavily (nearly verbatim) from Mandl and Shaw's Quantum Field 



12j. Starting with Maxwell's microscopic equations, 
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write the antisymmetric field tensor 

E x /c E y /c E z /c 
-E x /c B z -B y 
-Ey/c -B z B x 
-E z /c B y -B x 

where v £ {0,1,2,3} give the row and column, respectively and c = 1/^/io/Io is the 
speed of light in vacuum. The metric tensor is (700 = ~9jj — +1 f° r 3 {1)2,3} and 
9fj,u = otherwise, and = d/dx^. With the charge-current density s fl (x) = (cp(x), 3(x)), 
Maxwell's equations become 

d u F» v {x) = fios^x) , (4) 



d*F fMU (x)+d»F vX (x) + d u F 



. 



(5) 



Antisymmetry in F iiV gives d tl s il {x) = 0, conservation of the current coupled to the electro- 
magnetic field. In terms of the four-vector potential A^(x) = ($/c, A), the field tensor is 
expressed as 

F" v (x) = d v A»(x) - , (6) 

which is equivalent to 

d 

B = V x A , E = -V$ - — A . (7) 

dt v ; 

In terms of A^(x), Equations (j5J), known as the Bianchi identities, are satisfied identically, 
and Equations (jlj), known as "the" field equations, read 



d u d u {A»(x)) - d"{d v A"(x)) = Hos"{x) 



(8) 



which are Lorentz covariant and gauge invariant. Consistency upon quantization requires 
d l/ A u (x) = 0, which is recognized as the Lorentz condition, achievable from any given con- 
figuration by a gauge transformation. Incorporating "the plasma approximation" amounts 
simply to replacing the = signs with 7^ signs in Equations (T4|) and (jSJ). Thus, one sees 
that Gauss' law is not an isolated, arbitrary equation, being part of the unification of elec- 
trodynamics brought on by Maxwell [23], and that its violation amounts to decoupling 
magnetostatic fields from source currents as well as electrostatic fields from source charges, 
thus breaking the symmetries of the underlying gauge field theory, hence is incompatible 
with the laws of physics. 

Another way to look at the problem is to consider the manner in which students of plasma 
phys.cs are introduced to the concept (starting at the undergraduate level!) by following the 
arguments presented by Chen [1[ and Dougherty [16j (which in no way reflects a criticism 
of these particular authors, and in fact compliments them on such a clear elucidation of 
the topic at hand): Using continuity, take the divergence of Maxwell's curl equations to 
show that his divergence equations are "redundant" and may be dropped from the analysis. 
However, from the development above, we see that the natural, physical division of Maxwell's 
laws is not into the "curl" and "divergence" equations but rather into the homogeneous and 
inhomogeneous equations, and that dropping the electric field's divergence equation in favor 
of its curl equation amounts to the replacement of an information bearing field equation with 
one of the Bianchi identities, thus literally leaving useful information out of the subsequent 
analysis and hence losing touch with the real physics. At the classical level, simply rewriting 
the fields in terms of the four-potential obviates the question of redundant equations. The 
subsequent justification of the validity of "the plasma approximation" given by Chen in 
the context of the linear theory of ion acoustic waves actually only justifies the neglect of 
perturbations to the ion density down to scales on the order of the Debye length, not the 
neglect of Gauss' law, as its identical twin Poisson's equation is used throughout the analysis. 

Let us examine a recent appearance 24j and discussion 25|, [26] of the issue in the geophys- 
ical literature as it pertains to "large amplitude stationary structures in general plasmas" in 
a "static magnetic field" . Here, Verheest finds that demanding "strict charge neutrality [on 
a differential volume element, one must assume from the fluid formalism] has indicated that 
solitons and oscillitons cannot exist in electron-ion plasmas," rather than the suggestion by 
McKenzie, et al |27(], "that fully relativistic effects may prevent the formation of oscilliton- 
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like structures." From the discussion above, we see that these arguments are not unrelated, 
nor are they in conflict, as only by respecting Gauss' law does Maxwell's theory remain 
Lorentz covariant. The statement by McKenzie, et al, in their Comment, "that a question 
of using Poisson's equation versus the quasineutrality condition is not a new one [and] was 
discussed and clarified 40 years ago in the classical review paper by Braginskii |3]," is not 
accurate. The relevant discussion surrounding Equation (6.13) thereof is declaratory and 
unreferenced, hardly "clarified," yet has been taken as gospel by subsequent generations of 



plasma physicists in a variety of applications. Earlier treatments 4] focus on the argument 
that one cannot get more than one piece of information out of Poisson's equation, and as 
Equations ([S]) show, one does not, as the electrostatic potential and the supporting charge 
density are intimately related. 

Turning now to a popular application of "the plasma approximation" to the evaluation 
of the electrostatic field in a tokamak discharge at equilibrium within the prescripts of 
the neoclassical model 128 



(some geometric nomenclature is given in the caption 
of Figure 1, and these technical arguments serve merely to support, not to define, the 
main contention of this article), we note that using an impurity ion species' radial equation 
of motion to determine a non-vanishing, non-poloidally varying radial electrostatic field 
is at best an indirect measurement subject to the assumptions of the underlying model. 
(The direct measurement of the radial electric field accomplished by motional Stark effect 
diagnostics [3l|, [32] viewing the D a line emission of injected neutrals excited by collisions 
with plasma particles is biased towards observations along the horizontal midplane as a 
consequence of port placement hence does not necessarily reveal the true radial field profile in 
the (Z, R) plane; neither does it distinguish between static and dynamic electric fields.) The 
neoclassical prescription to use an equation of motion, rather than Gauss' law, to determine 
the radial electrostatic field, and then to turn around and use Poisson's equation to determine 
the charge density, is tautologically inconsistent, as the implication of V • E ^ p/eo must 
be that — V 2 $ 7^ p/eo- In his development of the fluid equation of motion from the first 
moment of the Vlasov equation 16], Elliot states that "the second integral [in Equation 
(2.11.10) thereof] is also zero, because E is not a function of" the particle velocities. Thus, 
neither is E a function of their average, the fluid velocity V. Rather, it should be considered 
as an input to the equation of motion, which is a differential equation to be solved for V. 
The defining relations for E are, naturally, Maxwell's, and the use of the fluid equation of 
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motion to determine E is inconsistent with its own derivation. 



giving its radial 



Nonetheless, the fluid equation of motion, Equation (5) in Reference 
component in the neoclassical model, must hold even when a velocity component vanishes, as 
does the radial component at equilibrium, leaving its interpretation in a bit of a conundrum: 
if the defining properties of the plasma state include the requirement of charge shielding, 
then the net charge on a differential volume element, which must be much larger than the 
Debye length yet smaller than any other scales of consideration for the averaging inherent 
in the definition of continuum quantities such as density, temperature, etc, to be formally 
valid, necessarily must vanish, as does any emitted electric flux. Taken at face value, the 
electrostatic potential, Figure 2, calculated from the radial field presented by Solomon, et 



al, in Figure 12(c) of Reference 28|, is not of harmonic form 33|, |34j, |35|, |36(, Figure 3, 
thus it implies a shifting of the electron density profile relative to the ions on the order of 
O(±10 13 /m 3 ), Figure 4, which does not sound like much when the absolute density is on the 
order of O(10 19 /m 3 ), yet produces a non-cancellation of terms which usually vanish in the 
quasineutral treatment, n e = J2ions z i n i f° r s P ec i es charge e s = z s e. The first inconsistency 
arises in the supposition of no poloidal variation to the radial field, as the radial equation 
of motion retains a poloidal dependence when the intrinsic variation of quantities induced 
by the geometry is considered [37|, given in the large aspect ratio limit, e = r / R 1, by 

B = B°/(l+ecos0) , V es = V e ° s /(l+ecos6) , V+, = V$ s (l + ecos9) , n s = n° s , (9) 

where A = (A) = § d6{l + ecos 6)A/2n is the flux-surface average of A. Using these values, 

Eriff) =P?/e s n° s + V° s B e -V e s B°/(l+Ecos6) 2 , (10) 

for species pressure p' s = dn s T s /dr, thus E r (9) ^ E®. (A further complication occurs if the 
species temperature T s = kT s retains a poloidal dependence.) Reconsidering that equation 
in light of a non-vanishing charge density, Yl s e s n s ^ 0, we sum over species to obtain 

E r e * n s = Yl & + esUs ( V <t> sBd - V 0s B $)\ = > ( U ) 

s s 

where the l.h.s. may be interpreted physically as proportional to 1 x 1 or 0x0 or neoclassically 
as 1 x 0. Noting that the r.h.s. is the radial component of the equilibrium equation Vp = 
J x B lets one identify it as zero, leading to inconsistency for a non- vanishing charge density. 
The alternative interpretation is to treat the plasma as a conducting fluid which is neutral on 



scales down to the Debye length and to reexamine the convective and viscous terms usually 
neglected, to wit, 



n s m s [(V a • V) V s } r +f ■ V-n s = n s e s (V 6s B^ - V^ s B e ) -p' s , (12) 

where we note that the convective term has a contribution ~ nmVg/r which survives incom- 
pressibility and the flux surface average. Thus, when one measures a species' net toroidal 
and poloidal velocities (in addition to the usual density, temperature, and magnetic field 
measurements), one might actually be determining that species' net radial viscous force. 

From our demonstration of the incompatibility of "the plasma approximation" with the 
laws of physics as expressed in the language of field theory, as well as our identification of 
some of the inconsistencies that arise from its implementation, what are we to conclude? 
Unfortunately and with some regret, we must conclude that much of the last half-century's 
worth of analysis on condensed matter in the plasma state within the fluid formalism of 
magnetohydrodynamics is gravely in error, being based upon beliefs and equations incom- 



patible with the laws of Nature. By neglecting Gauss' law, one cuts onesel 
powerful, and beautiful mathematics of field theory. Leslie C. Woods 38|, 



: off from the rich, 
39] was correct in 



his assessment of the state of fusion plasma research and the need to reevaluate some of the 
established beliefs long-cherished in the field. Perhaps with a more reactionary approach to 
plasma physics and fusion engineering, one which respects all of Maxwell's laws, the 21st 
century will see the fulfillment of "the promise of fusion energy 
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FIG. 1: Cylindrical coordinates (Z,R,<p) relate to tokamak coordinates (r,6,<j)) via Z = — rsin# 
and R = Rq + rcos9 in the concentric circular approximation, where r, 9, and 4> give the radial, 
poloidal, and toroidal directions, respectively. The outermost minor radius of the confined plasma, 
given in meters by a, defines the normalized minor radius p = r/a, and R a is its centroid. The 
magnetic field B and current density J lie in isobaric surfaces given by Vp = J x B for a stationary 
equilibrium, defining the "flux surface" at radius r. In general, the nested flux surfaces are neither 
circular nor concentric. 



Neoclassical Electrostatic Potential and its Derivative 
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FIG. 2: Neoclassical electrostatic potential <!> and its derivative E r = —d<&/dr vs normalized minor 
radius p. The potential is evaluated numerically from the radial field given in Figure 12(c) of 
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Neoclassical Electrostatic Potential 
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FIG. 3: (Color online.) Neoclassical electrostatic potential <3? in the (Z,R) plane, assuming no 
poloidal dependence, painted onto a generic equilibrium. Note that the potential is not of harmonic 
form, as the value at the center is not equal to the average of the values around the boundary, 
hence requires a supporting charge density. 

Electron Density Offset 
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FIG. 4: Flux-surface averaged electron density offset given by Poisson's equation required to sup- 
port the neoclassical electrostatic potential in units of 10 13 /m 3 . 
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